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Abstract. The phenomenon of shadowing of nuclear structure functions at small 
values of Bjorken-x is analyzed. First, multiple scattering is discussed as the underlying 
ly^ I physical mechanism. In this context three different but related approaches are 

presented: Glauber-like rescatterings, Gribov inelastic shadowing and ideas based on 
^kJ ■ high-density Quantum Chromodynamics. Next, different parametrizations of nuclear 

QQ I partonic distributions based on fit analysis to existing data combined with Dokshitzer- 

^^ ' Gribov-Lipatov-Altarelli-Parisi evolution, are reviewed. Finally, a comparison of the 

^^ , different approaches is shown, and a few phenomenological consequences of nuclear 

f-^ ' shadowing in high-energy nuclear collisions are presented. 

o 

:a 

> 

X 



Nuclear shadowing 2 

1. Introduction 

The fact that nuclear structure functions in nuclei are different from the superposition 
of those of their constituents nucleons is a well known phenomenon since the early 
seventies, see references in the reviews 012]. For example, for F2 the nuclear ratio is 
defined as the nuclear structure function per nucleon divided by the nucleon structure 
function, 






Rf2\^^Q I ~ A i7nuclcon/„ r\2\ ' \^) 



Here|, A is the nuclear mass number (number of nucleons in the nucleus). The variables 
X and Q^ are defined as usually in leptoproduction or deep inelastic scattering (DIS) 
experiments: in the scattering of a lepton with four-momentum /c on a nucleus with four- 
momentum Ap mediated by photon exchange (the dominant process at Q^ <C m'^o,m'^ 
where most nuclear data exist), 

lik) + A{Ap) ^ l{k') + X{Ap'), 

q = k-k\ W' = {q + p)\ x=^ = — ^^ , (2) 

see Fig. ^ The variable x has the meaning of the momentum fraction of the nucleon in 
the nucleus carried by the parton with which the photon has interacted. Q^ = —q^ > 
represents the squared inverse resolution of the photon as a probe of the nuclear content. 
And W"^ is the center-of-mass-system energy of the virtual photon-nucleon collision 
(lepton masses have been neglected and mnucieon is the nucleon mass), see e.g. [S] for full 
explanations. The nucleon structure function is usually defined through measurements 
on deuterium, ^°"'=i«^°'i = ^deutermm^2^ assuming nuclear effects in deuterium to be 
neghgible. 

The behaviour of Rp^{x,Q'^) as a function of x for a given fixed Q^ is shown 
schematically in Fig. El It can be divided into four regions§: 

• Rp^ > 1 for X > 0.8: the Fermi motion region. 

• i?4 < 1 for 0.25 ^ 0.3 < X < 0.8: the EMC region (EMC stands for European 
Muon Collaboration). 

• Rp^ > 1 for 0.1 < X < 0.25 ^ 0.3: the antishadowing region. 

• Rp^ < 1 for X < 0.1: the shadowing region. 

This review will be focused in the small x region i.e. that of shadowing, see PQ |2] 
for discussions on the other regions||. The most recent experimental data [U El El 13 IHl E! 

I Sometimes the ratio of nuclear ratios is used e.g. R{A/B) — Rp / Rp . 

§ Note that the deviation of the nuclear i^2-i'atios from one in all four regions of a;, is sometimes referred 
to as the EMC effect. I use this notation only for the depletion observed for 0.25 -;- 0.3 ^ a; < 0.8. 

II The region of Fermi motion is explained by the Fermi motion of the nucleons. For the EMC region 
there exist several explanations: nuclear binding, pion exchange, a change in the nucleon radius,. . . The 
antishadowing region is usually discussed as coming from the application of sum rules for momentum, 
baryon number,. . . 



Nuclear shadowing 




Figure 1. Diagram of leptoproduction on a nucleus through virtual photon exchange. 
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Figure 2. Schematic behaviour of Rp [x, Q^) as a function of x for a given fixed Q^ 



(see [H El Uni CU 1121 fo'^ previous experimental results), confined to a limited region of 
not very low x and small or moderate Q^ (and with a strong kinematical correlation 
between small x and small Q^, see Fig. O)), indicate that: i) shadowing increases with 
decreasing x, though at the smallest available values of x the behaviour is compatible 
with either a saturation or a mild decrease jH]; ii) shadowing increases with the mass 
number of the nucleus P; and iii) shadowing decreases with increasing Q^ jjj. On 
the other hand, the existing experimental data do not allow a determination of the 
dependence of shadowing on the centrality of the collision. 

In the region of small x, partonic distributions are dominated by sea quarks and 
gluons. Thus isospin effects, partially corrected in practice by the use of deuterium as 
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Figure 3. Plot on the left: Kinematical range in the x-Q^ plane probed in nuclear 
DIS gllSlElEIIHliniliniiniini and Drell-Yan 13 processes, and in d-Au at forward 
rapidities |14l I15| at RHIC. [Figure taken from |1^.] Plot on the right: The average 
values of x and Q^ of the DIS data from the New Muon Collaboration @1 [3 IHl C] 
(triangles) and E665 [HI E] (diamonds) in l-A^ and of X2 and AP of the Drell-Yan 
dilepton data \VA\ (squares) in p-^. The heavy quark mass scales are shown by 
the horizontal dashed lines. Those lines labeled saturation indicate the estimated 
saturation scale in proton and Pb. The different bands and lines show the values of x 
and Q^ which are or will be probed in Drell-Yan or heavy flavour production at SPS, 
RHIC and LHC, for rapidities different from central ones when indicated. [Figure 
taken from ^\-] See also the text in Subsection 12.31 and in Section El 



reference and of isoscalar nuclei, are negligible and will not be discussed in the following. 
In most approaches, the origin of the depletion of the nuclear ratios in this region is 
related with the hadronic behaviour of the virtual photon ^H] ■ This resolved hadronic 
component of the photon wave function at high collision energies - equivalent to small 
values of x, see (j21) - and at relatively low values of Q^, will interact several times with 
the different nucleons in the nucleus i.e. will experience multiple scattering. As I will 
discuss in the next Section, this results in a reduction of the corresponding cross sections 
- shadowing, related to the structure functions through 

with ttEM the fine structure constant. Thus, the phenomenon of multiple scattering is 
sometimes referred to as shadowing corrections. 

The importance of the phenomenon of nuclear shadowing is twofold: First, 
on the theoretical side it offers an experimentally accessible testing ground for our 
understanding of Quantum Chromodynamics (QCD) in the high-energy regime |19j . 
Multiple scattering is unavoidable in a quantum field theory as a consequence of such 
a basic requirement of the theory as unitarity. The nuclear size gives the possibility 



Nuclear shadowing 5 

to control the amount of multiple scattering at given values of momentum fraction 
X and scale Q. Besides, by varying the scale and the energy of the collision the 
interplay between the soft non-perturbative and the hard perturbative regimes can be 
addressed. Second, experimental studies on high-energy nuclear colhsions like those 
at the Relativistic Heavy Ion Collider (RHIC) [201 EB E21 ES] at the Brookhaven 
National Laboratory (BNL) are presently carried out. New possibilities like the Large 
Hadron CoUider (LHC) [HI at CERN or the Electron-Ion Collider (EIC) [24J under 
consideration, will become available in the future. They test the behaviour of parton 
densities inside nuclei at larger energies/smaller x than those presently available in 
fixed target studies like those at the Super Proton Synchrotron (SPS) at CERN. These 
new experimental data offer the possibility to further constrain our knowledge on 
the behaviour of nuclear cross sections and structure functions, both for observables 
characterized by a large scale for which standard perturbation theory can be applied, and 
for those with intermediate and small scales where new methods have been developed. 
The plan of the review is the following: In Section El models based on multiple 
scattering will be reviewed. In Section 01 those models which do not try to address the 
origin of nuclear shadowing but rather to study its evolution through the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations J2S1 EHl EI] , will be discussed. A 
comparison of the different models will be shown in Section HI Next, some consequences 
on high-energy nuclear collisions will be presented in Section Finally, in Section IHl 
some conclusions will be drawn. I have mainly focused on the most recent approaches 
roughly starting from the early nineties, please have a look at P for an extensive list of 
references on earlier theoretical and experimental work. 

2. Models based on multiple scattering 

As commented in the Introduction, the usual explanation for the origin of shadowing is 

multiple scattering [2HiEsiisniisiiisaissiisiiinsiinniinziisHiEniiiniiiiiii2iiisiiiaii3 

linillIlllHlliniEniEIlEaE3EllESI- WhUe the basis of the explanation is common, 
phenomenological details of its application vary from model to model. For example, 
the hadronic component of the virtual photon may be given a partonic structure like 
in the dipole model jSEl, see the next Subsection, or modeled [2HI EH HH IH21 IHB] as 
a superposition of hadronic states with the photon quantum numbers - vector meson 
dominance, or some combination of both approaches e.g. [3^. Besides, what is seen 
as multiple scattering in the rest frame of the nucleus corresponds to recombination in 
the infinite momentum frame jlHl EOl- Just to mention a few differences between the 
results of the models, the behaviour of models jHEl IHZI IHH] is dominated by hadronic 
configurations of large size, so the results turns out to be basically Q^-independent. On 
the other hand, models which consider an expansion in power-suppressed corrections 
in l/Q, either a fixed number of terms [KlU IKT| E21 or some re-summation [Sl|, show a 
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clear Q^ dependence^. Some models rely on an eikonal approximation jSHlllIlElj, see 
below, and are unable to reproduce the return of F2 nuclear ratios to 1 at x ~ 0.1, while 
others |12l ESI HH HEl UHl HH HHI include effects of finite coherence length, see below, 
and are able to reproduce such a behaviour. 

In this Section I will start by working out a little exercise which shows how multiple 
scattering leads to shadowing. This exercise should also clarify the origin of coherence 
effects. Then, in the Subsections models based on Glauber-like rescatterings, on Gribov 
inelastic shadowing, and finally the ideas based on high-energy QCD ^H], will be 
reviewed. 

For the exercise I consider the contribution coming from one and two scatterings, 
to the high-energy cross section of a massless scalar particle on a nucleus with mass 
number A. The scattering centers plus the interaction vertex are represented by the 
projectile-nucleon forward scattering amphtude times the nuclear density (see Fig. Hj). 
This example follows the spirit of the Glauber-Gribov theory [211 EH EH]; technical 
details can be found in Section 3.1 and Appendix A in [HHI for the case of scalar QCD at 
high energy. I will use light-cone coordinates a± = ao±az, a = (oq, a^, flz) = {0.+, 0.-, ^t) 



'J' 



(^^[(271)"* (k^+ie)]"" 
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Figure 4. One- (left) and two- (right) scattering diagrams, with the corresponding 
Feynman rules written on them. 



with ax = {ttx, tty) the two-dimensional transverse vector, assume dominance of the -I— 
components for the projectile, and define q = p' — p. I will employ the optical theorem 
for purely imaginary amplitudes, it{q = 0) = itforw = — o" for projectile-nucleon and 
iTn{q = 0) = — cr^ for the n-scattering contribution for projectile-nucleus collisions. 
Then the amplitude with one scattering (Fig. |3] left) reads: 

c{p+,p'^)iTi{q) = itior^ A(p+ +p'+) / d^xpA(a;+,XT)e'"-('''-P) 



% The experimental evidence of a Q^ dependence of the nuclear _F2-ratios comes from 0, thus being 
subsequent to many models e.g. the analysis in |56j : previous data did not show any clear Q^- 
dependence. 
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= ztforw cip+,p'^)A J d'xTTAixT)e-'''^<PT-P^\ (4) 



Pa{x+,xt) is the nuclear density normalized to 1, 



+ 00 



Ta{xt) = / dx+pA{x+,XT) (5) 



the nuclear profile, \xt\ = b the impact parameter and c{p+,p'_^) = {2n)2p^6{p'_^ ~ P+) 
a normalization factor. For the forward scattering case q = 0, ^ gives 

a\ = Aa (6) 

as expected. The two-scattering contribution (Fig. 0] right) reads 

c{p+,p'_,)z%{q) = zA{A - l)(ztforw)' I 0-,d'x^d'x2 e"-^'^-^) 

Xe ^^^ > — -— Pa{Xi+,Xit)Pa{X2+,X2t) 

= c(p+,p'+)A(A-l)(ztforw)' (7) 

X / f^dx,^dx2+d^X,Td\2T e-^4(-2+--i+)/(2p+) 

J (27r)^ 

with 9{x) the step function. The second equality follows from doing the dk- integral 
closing the integration contour with an infinite semicircle in the lower half-plane and 
simplifying the remaining (5-functions. 

Coherence effects are contained in the factor exp [— 'i/i;|,(x2+ — Xi+)/(2p+)] which 
can be written as exp [—i{x2+ — Xi^)/l^, with l^ = 2p+/k^ the coherence length. In 
the low energy limit p+ ^ this factor leads to iT2{q) —>■ (the same happens for 
contributions with more than two scatterings). Then all rescattering corrections vanish, 
so the total cross section results equal to the superposition of single scatterings (jHI) - 
this limit is thus called the incoherent limit. This would imply a nuclear ratio equal to 
1, so shadowing vanishes for low energies i.e. large values of x. On the other hand, in 
the large energy, totally coherent limit j9+ ^ oo, exp [—i{x2+ — xi+)//c] — ^ 1 and (jZj) 
gives 

tT2iq) = ^i^^{it,,,„Y I d'xT e— (p^-^-)tI(xt), (8) 

which in the forward case gives 

a\ = - ^^\~^^ J d'xT [TAixryr. (9) 

This correction turns out to be negative so the cross section is smaller than the 
superposition of independent collisions of the projectile with every nucleon in the 
nucleus, and the nuclear ratio results lower than one. In this way, multiple scattering 



Nuclear shadowing 8 

offers an explanation for shadowing"'". Furthermore, from this example it becomes 
evident that shadowing increases with increasing a and mass number A (and, if the 
integration over xt is not performed, increasing centrality equivalent to decreasing 
xt)- The cross section a increases with increasing energy (equivalent to decreasing 
x) or decreasing Q^ (equivalent to increasing size of the hadronic component of the 
virtual photon, see the next Subsection). Most of these features, as indicated in the 
Introduction, are seen in the experimental data. The differences between the models 
mentioned at the beginning of this Section come both from the modeling of a and to 
the way in which multiple scattering is considered. Both aspects are the subject of the 
following Subsections. 

In the coherent limit the hadronic fluctuation of the virtual photon interacts with 
the target as a whole. The kinematical region where this happens can be discussed 
using a simple argument: In the laboratory frame where the nuclear target is at rest, the 
lifetime of the hadronic fluctuation of a 7* with squared virtuality Q^ can be estimated 
using the uncertainty principle and the Lorentz dilation in this frame, 

1 i^lab W^ 1 ..„. 

T ^^ — X ~ ~ (10) 

^/ ^/ -^ ''^nucleon ^5 ■^''''nucleon''' 

where I have used (0) for 2mnucieon-Ezafe — W^ ^ Q^. This lifetime increases with 
decreasing x or increasing energy. For the hadronic fluctuation to interact with the 
nuclear target as a whole, the lifetime has to be greater than the nuclear radius, 
r > Ra- This implies x < l/(2mnucieon-Ryi)- Using typical values Ra ~ A^^^ fm, 
we get X < 0.1A~^^^ which roughly coincides with the experimentally measured values 
of X for which the transition from antishadowing to shadowing takes place. 

2.1. Glauber-like rescatterings 

Some models try to address the origin of nuclear shadowing through the Glauber-Gribov 
formalism in the totally coherent limit [311 El lEH] • A proper treatment of coherence 
effects requires the consideration of the mass spectrum of intermediate fluctuations and 
will be discussed in detail in the next Subsection, though some models jlHI include 
coherence in an effective way. 

The Glauber-Gribov theory [SZl IHHl EHj considers the multiple scattering of the 
hadronic component of the virtual photon with a nucleus made of nucleons whose binding 
energy is neglected. This hadronic component keeps a fixed size during the scattering 
process - the eikonal approximation, and is usually limited to its lowest-order Fock state, 
a qq pair - the so-called dipole model jSHlinSinni- Then the total dipole-nucleus cross 
section reads 



o'dipoic- A{x,r) = d% 2 



exp ( --ATA{b)adipolc-nnc\eon{x,r) 



(11) 



"•" Multiple scattering plays a key role in many physical processes. For example, coherence effects in 
multiple scattering are widely discussed in the context of medium-induced radiation both in Quantum 
Electrodynamics and in QCD (see e.g. |f)lj and |()2| respectively, and references therein), or in heavy 
flavour production on nuclear targets |63[ IM| . 
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with h the impact parameter of the center of the dipole relative to the center of the 
nucleus and r the size of the dipole. This cross section is then related to the nuclear F2 
through (jHl) and 

(y^,_A{x,Q'^) = / d'^rp{r,Q^) adipoic- A{x,r), (12) 

where p(r, Q"^) are the distributions of colour dipoles of size r created by sphtting of 
the incident photon into a qq pair ISEl ESI EEI • These distributions provide a definite 
relation between increasing Q and decreasing r (or vice versa). 

The differences between realizations of this approach come mainly from the model 
used for the dipole-nucleon cross section (equivalent to the dipole-proton cross section 
at the small values of x where this approach is applicable). For example, in |1T] 
a parametrization based on a saturation model [HZ] is used, while in jSH] a form 
based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pomeron j68[ lU^ is employed, see 
Subsection 12.31 On the other hand, in [lOj forms based on the double-leading-log (DLL) 
approximation to the DGLAP evolution equations are taken, see jTUI for a discussion 
of the relation of scales in the dipole model and DGLAP. All of them give a reasonable 
description of the data on nuclear shadowing for x < 0.01, see e.g. Fig. although the 
effective introduction of coherence effects in [IHj allows to describe the whole shadowing 
region. The Q^-dependence of nuclear ratios [Z] is also reproduced jH]. It must be 
stressed that once the dipole-nucleon cross section is fixed, the extension to the nuclear 
case is essentially parameter- free, making the agreement with the nuclear experimental 
data more remarkable. 

Apart from their intrinsic interest (see e.g. an application to exclusive vector meson 
photoproduction in |7Tj), these models have also been used as initial conditions at not 
very small x ~ 0.01, for evolution towards smaller values of x in the framework of high- 
density QCD, see Subsection 12. HI Also there I will discuss the issues of the saturation 
scale which can be extracted in this framework. 

2.2. Gribov inelastic shadowing 

In the classical Glauber model |2ZI subsequent interactions of the projectile with 
nucleons in the nucleus occur, and the intermediate states of the projectile are the 
same as the initial one i.e. elastic. In the relativistic Gribov theory [HSl 1^ - 
subsequent interactions are suppressed at high energies and the collision proceeds 
through simultaneous interactions of the projectile with the nucleons in the nucleus. The 
intermediate states are no longer the same as the initial state and are called inelastic. 
The use of Reggeon calculus [2! and the Abramovsky-Gribov-Kancheli (AGK) cutting 
rules [Zni (see an updated discussion in [Hj) allow to write a relation between the cross 
section for diffractive dissociation of the projectile and the two-scattering contribution 
to the projectile-target cross section, see Fig. IHl 

The first correction to the non-additivity of cross sections comes from the second- 
order rescattering a\. In Fig. El diffractive DIS is shown both in the infinite momentum 
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Figure 5. Plots on the left: Nuclear size dependence of the i^2-ratios in the Glauber 
model in ^Ij for two fixed values of x, compared with experimental data ^: (filled 
points). Plots on the right: Q^-dependence of the nuclear i^2-ratios for two fixed 
values of x, compared with experimental data |Zj (filled points). Open triangles and 
circles in the plots on the left, and solid and dashed lines in the plots on the right, 
correspond to different models for the dipole-nucleon cross section. In the experimental 
points, inner error bars correspond to statistical errors, and the outer ones to statistical 
plus systematic errors added in quadrature, [Figures taken from 41 ,] 






Figure 6. Diagrams relating diffraction with the two-scattering contribution to the 
total cross section. 



frame and in the rest frame of the nucleon. From Fig. IHlit becomes clear that the square 
of such contribution is equivalent to a double exchange with a cut between the exchanged 
amplitudes, a so-called diffractive cut. To compute the first contribution to nuclear 
shadowing a\ which comes from these two exchanges, one needs its total contribution 
to the 7*-nucleon cross section. This contribution arises from cutting the two-exchange 
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amplitude in all possible ways: between the amplitudes and the amplitudes themselves 
in all possible manners. For purely imaginary amplitudes, it can be shown [7^ 175] that 
this total contribution is identical to minus the contribution from the diffractive cut. 
Thus diffractive DIS is directly related to the first contribution to nuclear shadowing. 
The final expression reads 



0". 



d%Tl{h) 

Jm'2 . 

rmn 



dM' ''-" 



dM^dt 



Fj^{tmin), (13) 

i=0 

.. , ...... - 4m2 and M^„^ = 

Q^ {xpmax/x — 1) with xpmax ~ 0.1, See jlH] . The usual variables for diffractive DIS: 
Q^, X, M^ and t, or xp = x/(3, (3 = Q^/{Q'^ + M^), are shown in Fig.[7| da'^,_^/dM^dt is 



with M^ the mass of the diffractively produced system, M^ 



. lepton ^ ^^^' 
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Figure 7. Diffractive dissociation in the infinite momentum frame (IMF), with the 
corresponding kinematical variables, and in the rest frame of the hadron. 



the differential cross section for diffractive dissociation of the virtual photon. Coherence 
effects are taken into account through 



FAitmin) — d b JoiW—tmin)TA 



(b), 



(14) 



with tr 



-m 



nucleon P 



This function is equal to 1 at x ^ and decreases with 
increasing x due to the loss of coherence for x > (2mnucieon-Ryi)~^, see f|Tnjl . Details can 
be found in the corresponding references jl21 113 EH HSl EHl ITTIIIH] . 

The differences between available realizations come from the consideration of real 
parts in the pomeron amplitude |10l EZ] , or from the model or parametrization used for 
the diffractive cross section: phenomenological models and parametrizations jlHl HEl EZl 
I3H1 which reproduce the existing experimental data*, or a model which considers the 
evolution of a diffractive qq state in |1^ 1^ Elj (see the method in [7^). Differences 

* The use of such parametrizations can be justified at large enough Q^ by the factorization theorem 
for diffractive hard scattering |75| . 
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also come from the extension of the models to include higher order rescatterings. Such 
extensions are mo del- dependent. Some models consider that all intermediate states in 
the subsequent rescatterings have the same structure [13 EEl IHl HE] , in the form of a 
Schwimmer ^^ 

^Sch _^ f r2r ^TAJb) 

or an eikonal unitarized cross section 

.^k _ f .2. A{1- exp [-2(A - l)T^(&)/(x, Q')]} 

(with f{x, Q"^) defined to get consistency with (fT^ when these equations are expanded 
to second order). Other models take into account the possibility of different intermediate 
states Ha 113 111- 

Several comments are in order. First, once the model or parametrization for the 
diffractive cross section is provided, the extension to the nuclear case is parameter-free 
- except for the modeling of higher order rescatterings (|T5|l and ()16|1 . These models, 
together with those in Subsection \2.\\ provide an impact parameter dependence of 
nuclear shadowing. In this way, this approach and the one presented in the previous 
Subsection offer a link between the nucleon and nuclear cases. Second, these models 
can be used as initial conditions for DGLAP evolution as done in jlHl IH] , see the next 
Section. Third, both models for the diffractive cross section and their extension to the 
nuclear case |^l^ do not correspond to any definite order in a power expansion in 1/Q 
but perform a re-summation of all powersjl. Finally, in these models the comparison with 
experimental data, when such comparison is available, turns out to be reasonable, see 
Fig. ISl- although the Q^-dependence of the nuclear ratios results too smooth compared 
to data [7 , which apparently indicates the need of additional DGLAP evolution. 

2.3. High- density QCD 

High-density QCD - the domain of large gluon densities - has become a very fashionable 
subject in the last fifteen years. It deals with the behaviour of QCD at very large 
energies. Regarding the contents of this review, it offers a definite theoretical framework 
to compute shadowing corrections although the high-energy approximations involved 
make its applicability to the present experimental situation a subject of intense debate. 
The literature on this matter is vast and I will not cite but a few papers, referring the 
reader to the contributions in ^^ and to the recent reviews [ZHl IZHl lEDl EIj where all 
the relevant references can be found. 

jj In 0B1 07] the discrepancy between the data and the results of the model when evolved through 
DGLAP to smaller values of Q^ from that, Qq = 4 GeV^ which is taken as initial value for evolution 
and where the parametrization of data is used, is considered as evidence of the existence of large power- 
suppressed contributions. On the other hand, in other models for diffractive data 021 EHl the presence 
of strong Q^-dependent terms is not required to describe the nuclear data at low Q^ . 
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Figure 8. x-dependence of the nuclear F2-ratios in tiie model in j48j for different 
nuclei, compared with experimental data |S1 IHI (circles). Solid lines correspond to 
Schwimmer (I15|) and dashed lines to eikonal (|16|l unitarization. Error bars in the 
experimental points follow the same convention as in Fig.[Sl In the ratios Xe/D, filled 
circles correspond to the analysis with hadron requirement and open circles to that 
with electromagnetic cuts, see the experimental paper |H] for more details. Both the 
experimental results and the theoretical ones, joined by lines, correspond to different 
average Q^ for every different value of x. [Figure taken from |48j.] 



In high-density QCD the small x partons (slow gluons) are treated classically due 
to the high occupation numberff oc l/a^. This number is as high as it can be - thus 
this field is often referred to as saturation physics. The source term for the classical 
equations of motion comes from the fast partons e.g. valence quarks with large x (see 

ft For example, in the BFKL framework |68[ I69| the gluon density xg in the hadron is expected to 
increase with decreasing x, ex x"^'^^"'. The exponent in this power takes a value ^ —0.5 for the 
strong coupling constant as ^ 0.2. DIS proton data show an increase ~ x~^^^ for small x, although no 
conclusive evidence of BFKL dynamics has been extracted from such behaviour. 
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Fig- El)- The gluon density per unit impact parameter and transverse momentum of the 
gluon, the so-called unintegrated gluon density at fixed impact parameter, computed 
in this way for an ultra-relativistic large nucleus - the McLerran-Venugopalan (MV) 
model E21 EHl El - reads 

_ d{xgA) 1 [ d'^XT 



dNf 



dyd^hd^kT d^hd^kx 



oc 



a.o 



-ixT-kx 



Jbn 



(^l _ e-4Q?/4J . (17) 



Q^g oc ATA{h)xg^^c\con is the squared saturation momentum or saturation scale, which 
corresponds to the typical gluon transverse momentum and to the scale at which the 
exponential in ()17|) starts to give large corrections. This saturation scale increases with 
increasing nuclear size and increasing energy or decreasing x. So it is plausible that 
at some given high energy, a^ becomes small enough for perturbative methods to be 
applied 
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classical 



Figure 9. Diagram showing the separation between fast and slow partons, and the 
contribution from radiated gluons. 



With increasing energy, radiation processes start to contribute, see Fig. IHl 
Additional gluons are radiated from the source in a kinematical region intermediate 
between the fast and slow ones, and are absorbed in a redefinition of the source. 
Mathematically this procedure results in a renormalization group equation for the 
distribution of colour sources in the hadron. Under several simplifications, this 
renormalization group equation gives a single closed non-linear equation for the dipole- 
hadron scattering amplitude^, the Balitsky-Kovchegov (BK) equation jSHUHnj. The use 
of the dipole model, see Subsection 12. ![ provides its link with the nuclear structure 
functions. Note that in the MV model, the number of partons is not modified but they 
are only redistributed in transverse momentum, while non-linear BK evolution does 
diminish the number of gluons. 

From the explicit form of the MV model ()17|) . it is obvious that it corresponds 
to a Glauber-like re-summation of rescatterings in the totally coherent, high-energy 

I See Hll(l . (|12|l : its relation with the unintegrated gluon density comes through the Bessel-Fourier 
transform defined in the right-hand side of H17|l . 
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limit. Indeed, the MV model can be used as initial condition at some not too small 
X ~ 0.01 for evolution towards smaller x through the BK equation. Other initial 
conditions have been essayed in the literature. But it is a noticeable property of the 
BK equation that all initial conditions result in a universal form for the solution of 
the equation after large enough evolution [HZl ESj- Such scaling implies, through the 
dipole model, that virtual photon- nucleon cross sections are not a function of x, Q^ 
and nuclear size separately, but only of Q^ IQ\ where all dependence on x and nuclear 
size is included in Qg. This scaling has been found in leptoproduction data on nucleon 
and nuclear targets jHH IHDl HI] • Nevertheless the situation is not yet clear: numerical 
solutions [021 EB] of the BK equation show that the asymptotic scaling behaviour appears 
at rapidities larger than those presently available, which are still dominated by the initial 
conditions. Besides the effects of a running coupling, not included in the derivation of 
the equation, are large. 

Predictions [^3 113 IHE] for nuclear structure functions of large nuclei at very small 
X have been computed in this framework. Some of them are presented in Fig. E3 
As the main results in this approach, the nuclear structure function at very small x 




- 




U" 


J} 


X 
X 


Z<0' 






1 


-^Q 



'0 10 10 10 10 



-0 10 10 10 10 -0 -0" 10 



Figure 10. Plot on the left: F2 for Pb at small x for different Q^ computed [SS] 
through the BK equation and the dipole model. Plot on the right: nuclear dependence 
a, F/(x, Q^) = A"(^''5^)i^^=i(a;, Q^) versus x for different Q^. [Figures taken 
from P3].] 



increases as In^ x with decreasing x. It also shows a Q^-dependence weaker than oc Q^ 
and a very strong nuclear dependence, with a(x, Q^) (F^(x, Q^) = A"^"^'*^ )F^^^(x, Q^)) 
getting smaller than the geometrical 2/3 factor for very small values of x §. Different 

§ Both this very small value of a and one of the powers of In x have their origin in the contribution of 
the dilute nuclear surface where the non-linear term in BK evolution, or saturation effects in general, 
are negligible. 



Nuclear shadowing 16 

realizations |M1 [HHl IHE] differ in the initial conditions, in the consideration of impact 
parameter, and in the treatment of large-size dipoles and of the evolution for not very 
small X. They turn out to give results which may vary as much as a factor 10 for 
X ~ 10~^. Predictions for heavy flavour production also exist jHSHlZj- 

The saturation scale computed within BK evolution behaves like Qs oc x~°'°% with 
(i = 4 -^ 5. Its dependence on the nuclear size is not yet fully determined; in the 
most widely employed approximation valid for a very large nucleus, the A-dependence 
follows that of the initial conditions, usually oc A^^'^. Besides, running coupling effects 
modify both dependencies dramatically [HSI- The saturation scale can also be studied 
within phenomenological approaches [HI EH lEHl IHIIl IHI]- For example, a value for 
the saturation scale can be obtained from Glauber approaches flTT|) as the value of Q^ 
for which the effect of the exponential factor in this equation becomes sizable (other 
geometrical criteria have also been essayed, like percolation [98j). Values extracted from 
this kind of studies are Ql ~ A'^(a;/0.01)-°-3 GeV^, with 5 > 1/3. 

Finally, other approach to the problem considers power-suppressed corrections || 
in 1/Q. Such power-suppressed contributions are enhanced by the nuclear size. The 
first power-suppressed correction to DGLAP evolution |HIH IKT| results in a non-linear 
equation. From the equality of the linear and non-linear terms, a value for the saturation 
scale can be extracted ^7] which results in rough agreement with the estimations 
previously discussed, see the solid black lines in Fig. El [Ej. More recently, such 
power-suppressed contributions have been re-summed |31] in the high-energy eikonal 
limit, resulting in a rescaling of the x variable whose results reasonably describe the 
experimental data, see Fig. ^2 Besides they are in agreement with available data on 
the nuclear effects on the longitudinal to transverse cross sections |innj . Also more 
phenomenological studies [SH] are in agreement with the experimental data. 



3. Models based on DGLAP evolution 

Another type of models do not try to address the origin of nuclear shadowing (or of 
modifications of parton densities in nuclei in general) but to study the Q^-evolution of 
nuclear ratios of parton densities, 

through the DGLAP evolution equations J2S1 I2H1 I2ZI , see also 0. From the very first 
attempts [M], several analysis have appeared HEl iZl EHl EH CUSl Ell CUSl CUE] ■ They 
try to perform for the nuclear case the same program developed for the nucleon: Nuclear 
ratios are parametrized at some value Qq ~ 1 ^ 2 GeV^ which is assumed large enough 
for perturbative DGLAP evolution to be applied reliably. These initial parametrizations 
for every parton density have to cover the full x range < x < 1. In the nuclear case, 

II The high-density QCD approach does not correspond to a fixed order in the power expansion but 
re-sums, in some hmit, aU power-suppressed contributions. 
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Figure 11. x-dependence of the ratios in the model in [SJ for different nuclei, 
compared with experimental data 0I|S1 121121 (filled points). The band corresponds to 
different choices of the scale of power corrections ^^ in |S31. Ai3_t = Data — Theory, 
and the open circles joined by dashed lines in these plots show the comparison to 
the approach in |99| . The (x, Q^) correlation of the experimental points is taken into 
account in the theoretical results shown here, as it was in those in Fig.|Sl [Figure taken 
from 15^.] 



the nuclear size appears as an additional variable. Then these initial conditions are 
evolved through the DGLAP equations towards larger values of Q^ and compared with 
experimental data. From this comparison the initial parametrizations are adjusted. 
Different approaches differ in several details, see |T7j: 

• The form of the parametrizations at the initial scale. For example, in [HHl ll()''?j 
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shadowing saturates for very small x^, contrary to |l(Jtij . Also the value of Ql 
varies e.g. from ~ 0.4 GeV^ ^^ to 2.25 GeV^ [111 [103 • The parametrizations 
for sea quarks and gluons in |l()4j do not show any EMC effect. Special mention 
has to be done to the approach of ^Hl ^\ where the initial gluon density is taken 
from diffractive nucleon data at Qg = 4 GeV^ as discussed in Subsection 12.21 and 
no attempt is made to modify it from the comparison with experimental data. 

• The use of different sets of experimental data. For example, Drell-Yan data |^ are 
used in jHH I1U3| llU5t I1U6J but not in |lU4j . These data give the main constraint 
to the valence and sea contributions in the antishadowing region in [HHl HO-'^j but 
the parametrizations in |105j do not show antishadowing for sea quarks. HERMES 
data |l()8j are used in |l()5j . Also the data on the Q^-dependence of nuclear ratios [Ij 
are included in [^ llU3t I1U5| I1U6J but not in jlU4j ; they give the main constrain on 
the gluon distribution at small and moderate x, see below. 

• The order of DGLAP evolution. The evolution is made at leading order (LO) 
in innmUS! and at next-to-leading (NLO) order in [HHl HHSl HUE] • This turns out 
to modify the Q^-dependence of nuclear ratios. 

• The treatment of isospin effects and the use of sum rules as additional constraints 
for evolution. For example, isospin symmetry of the nuclear ratios is assumed 
in |99 | llOHj but not in |l()4j . Momentum, charge and baryon number conservation 
are used in |1U4|I1U6J . but charge conservation is not used in j^llU3j . In practice, 
these differences result numerically small. 

• The different nucleon partons densities used in the analysis. In practice this choice 
is of little importance at the level of the nuclear ratios, as its effect appears in both 
the numerator and denominator in (fTH)) and cancels to a large extent. 

A comparison of different approaches can be found in Fig. El see also Section 0] 
The differences are noticeable, even more when one considers that all approaches have 
been designed to reproduce available experimental data (but see below the discussion 
on the Q^-dependence of nuclear ratios), see ^Zj for comments. Concerning shadowing, 
it turns out that the nuclear ratios for gluons are almost unconstrained for x < 0.02, 
and for sea quarks for x < 0.005. The stronger constraint on gluons in the region 
0.02 < X < 0.2 comes from the Q^- dependence of nuclear ratios which I will discuss 
now. 

The DGLAP equations establish a relation between the logarithmic Q^-evolution 
of the structure functions and the gluon distribution |109j . Such relation, valid at LO 
and small x, has been extended to the nuclear ratios |ll()j : 

dRJi^^Q') ^ lOa. xg{2x,Q') f ^ , _ ^ ,^ 

% The model |1U7) uses Glauber-like rescatterings to effectively include some high-density corrections 
to the parametrizations [991 lll)3| . These corrections lead to a shadowing which docs not saturate at 
small X, and which increases with increasing Q^ for small x < 10^"'. 
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Figure 12. Nuclear ratios for Ca versus x computed for different partons densities, for 
two different values of Q^. Results from |106| (solid), |S9| (dashed) and 104 (dotted 
lines) are compared. [Figure taken from |106 |.] 



In this way, Rf{2x,Q'^) > Rp,^{x,Q'^) implies a positive Q"^ slope, while Rf{2x,Q'^) < 
Rp^{x,Q^) gives a negative slope. The available data on the Q^-dependence of 
nuclear ratios |7j allow to constraint within the DGLAP evolution scheme, the relation 
between the nuclear gluon distribution and the nuclear ratio for F2. In Fig. ^J 
the result of Q^-evolution of the nuclear F2-ratio for Sn over C is shown |ll()j for 
different models [HHl I1U4J . and also for DGLAP-evolved parametrizations |11H I112J 
(these parametrizations were originally proposed as Q^-independent). While in the 
parametrization |112j the ratio for all flavours is equal at the initial scale Qq, and thus 
it gives a too small but positive slope, in the parametrization of jlllj the shadowing 
for gluons is much larger than that for sea and valence quarks so it results in a negative 
slope at the smallest x and Q^. Therefore, from the comparison with experimental 
data j2|, DGLAP analysis favours those sets in which gluons are less shadowed than 
quarks for x ~ 0.01. This is at variance with some approaches e.g. [JOlEl where gluons 
are more shadowed than quarks, see the next Section. As commented in Subsection I2.2| 
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the discrepancy between the data and the results of this model when evolved to smaller 
values of Q"^ from Ql = A Ge^ 
power-suppressed contributions^ 



values of Q^ from Qq = 4 GeV^ is considered as evidence of the existence of large 




Figure 13. Results of the models IHSl (solid) and |1()4| (dashed-dotted) , and of 
DGLAP evolution of initial conditions from |112| (dashed) and |111| (dotted lines), 
for the Q^-evolution of the ratio Sn/C at different values of x, compared with the 
experimental data [7]. [Figure taken from |11(J) .] 



While DGLAP approaches do not address the fundamental problem of the origin 
of shadowing, they are of great practical interest. They provide the parton densities 
required to compute cross sections for observables characterized by a hard scale for 
which coUinear factorization |113j can be applied, see Section El and [TTJ for discussions 
and further references. As a final comment, the centrality dependence of shadowing 
is not addressed in these models as the existing experimental data do not allow its 
determination, although some approaches e.g. |11H I114J provide an ansatz for such a 
dependence. 



4. Comparison among different models 

In Fig. El a comparison of the results of different models for R^ is shown. The results 
coincide within ~ 15 % in the region x ~ 0.01 where experimental data exists. But 
they strongly disagree for smaller values of x, the difference being almost a factor 2 

"*" DGLAP equations consider only leading power-suppressed contributions - thus the shadowing 
produced by DGLAP evolution is referred to as leading-twist shadowing. Nevertheless attempts have 
been done to include the first power-suppressed corrections |5()ll51j in this analysis of the Q^-dependence 
of nuclear ratios jlOlllllO] . But the result of such contributions is to make the logarithmic slope smaller, 
so they do not help to solve this problem. 
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at X ~ 10~^. In general, models ^Hl SH] based on Gribov inelastic shadowing give a 
larger shadowing than those based in Glauber-like rescatterings 1101 1^ . The model [02] 
based on high-density QCD gives less shadowing. Among the models based on DGLAP 
evolution, jOH I1U3J gives larger shadowing than |lU4j and to the one in |lU6j , see Fig. [T21 
for sea quarks which determine F2 for small x. But one should take into account that 
in DGLAP approaches the small x behaviour at Q^ close to Qq comes mainly from 
the assumptions for the initial conditions. The model |107j gives results for R^ at 
Q"^ = 2.25 GeV^ very close to those of [001 CHS! ^^ ^ — O-O^ where the considered high- 
density corrections are almost negligible, but sizably smaller, R^ ~ 0.43, at x = 10^^. 
These predictions could be checked in the EIC [21], where values of x ~ 3 • 10~^ at 
Q^ ~ 1 GeV^ should become accessible, or in ultra-peripheral proton-ion or ion-ion 
collisions (UPC) fTTBllrifi] . 
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Figure 14. F2-ratios for Pb versus x at fixed Q^ = 3 GeV^ from tiie models: Armesto 
et al. ("this worlc") |3H|, HKM [TDll, Sarcevic (SOj, Bartels ^, Frankfurt [4^ (at 
g2 ^ 4 GeV2), Armesto [Hj and EKS98 [HSllinSj- [Figure taken from gHj.] 



Now let me turn to the nuclear gluon density. In those approaches |4m HT] which 
rely on the dipole model, the gluon density is obtained from the unintegrated gluon 
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distribution through 



2 



.,,,(.,0^)^j%^krj,%§^^, (20) 

where A^ is some infrared cut-off, if required. This identification is only true for 
Q^ > Qsat W^ EHl CHI- In models based on Gribov inelastic shadowing US 11201, 
the nuclear gluon density is obtained from (fT ^ - (fT^ but using the gluon density in the 
pomeron instead of the diffractive structure function of the proton. 

In Fig. El a comparison of the results of different models for E^^ is shown. At 
variance to the case of -R^^, now there are large discrepancies also at a; ~ 0.01, see the 
discussion on the Q^-dependence of the nuclear ratios at the end of Section El This is 
due to the fact that the gluon density is only indirectly constrained by experimental 
data - F2 is mainly determined by the sea quarks at such values of x. The discrepancy 
between the different approaches is roughly a factor 2 at x ~ 0.01 and a factor 3 at 
X ~ 10~^. Among the DGLAP approaches the results of |lllj . in which gluon shadowing 
is fixed to reproduce the multiplicity in Au-Au collisions at RHIC (see the next Section), 
are clearly below those from [HT H llOHj and |l()4j *. as discussed at the end of Sectional 
Again, these DGLAP results are mainly determined by the initial parametrizations. The 
results of Glauber-like models are similar at x ~ 10~^ but those of jH] are smaller than 
the ones in jlOj at x ~ 0.01. Results from Gribov inelastic shadowing jlB] show a strong 
shadowing: see Fig. Upland the mentioned discussion in Sectional |120j gives E^^ ~ 0.7 
at X = 0.01 and ~ 0.2 at x = 10"^ for Q^ = 6.5 GeV^. The model ^\ gives resuhs for 
i?Pb ^^ Q2 ^ 2.25 GeV^ very close to those of ^ (103 at x = 0.01 - as it was the case 
for Rf^, but sizably smaller, i?^^ ~ 0.37, at x = 10"^ Finally, |l2j gives Rf" ~ 0.95 at 
X = 0.01 and ~ 0.75 at x = 10"^ for Q"^ = A GeV^. 

Contrary to F2 which is directly measurable, the gluon distribution is only indirectly 
constrained both at an EIC j21|, in proton- nucleus or nucleus-nucleus collisions at 
RHIC |2niinil21l2Sl and the LHC [HI, and in UPC [113 Eni. Some of these indirect 
constraints in proton-nucleus and nucleus-nucleus collisions will be discussed in the next 
Section. 

5. Shadowing and particle production in high-energy nuclear collisions 

In this Section I will discuss some consequences of the phenomenon of shadowing in high- 
energy nuclear collisions. I will devote most of the Section to proton(deuteron)-nucleus 
collisions, and only at the end I will briefly refer to nucleus-nucleus collisions. 

Particle production is related to the parton densities in the projectile and target 
through a factorization theorem, which in general takes the schematic form 

aA-B^CX OC (t)i/A{xA) ® (t)j/B{xB) ® CFij^c{xA, Xb, Qc)- (21) 

* This model gives results very close to those from |106| . see Fig. ^1 for gluons, although the latter 
shows an increase towards one at very small x due to the form of the initial parametrizations. 
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Figure 15. Ratios of gluon distribution functions for Pb versus x from different 
models at Q^ = 5 QeV^: HKM jTOl, Sarcevic gO], EKS98 [23110^1, Frankfurt gSI, 
Armesto [H] and new HIJING |lllj . RHIC and LHC point the ranges of a; = 



{Q/y/^)ey for processes with \y\ < 0.5, Q^ = 5 GeV^ at RHIC (^/s^^) 
GeV) and LHC (ys^) = 5.5 TeV) respectively. [Figure taken from XL-] 
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In this formula, integration over the relevant variable is implicit, a is a scattering 
matrix for partons i, j to give the produced parton/particle C. Qc is a scale (e.g. mass, 
transverse momentum,. . .) characteristic of C. And (pi/AixA) are the probabilities of 
finding parton i in hadron/nucleus A with momentum fraction xa- Depending on the 
type of factorization, collinear for Aqcd ^ Qc ^ a/s I113J or kx for Aqcd ^ Qc ^ 
\fs |12H II22I I123J . the matrix elements are different and the partons distributions are 
integrated or unintegrated respectively |119j . For the status of collinear factorization in 
collisions involving nuclei, see ^]\- The validity and exact form of kx factorization in 
proton-nucleus collisions - here proton means a hadron in which the parton density is 
small - is under discussion, see e.g. |124[ 1125't I126J and references thereinjj. For nucleus- 



tt In |127[ I125| . a non-linear fc^ factorization is proposed and the properties of the resulting nuclear 
gluon density are studied. 
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nucleus |129t 11301 H^lj , arguments exist which indicate that it is not vahd |132t I133J , 
but the size of the violations has not been fully quantified yet. Therefore, most of the 
discussions in the literature assume the validity of such factorizations in the nuclear 
case, which will be done in the following. 

In a 2 ^ 2 process the relation between rapidity y and transverse mass itit — 
^Jrn? + p^ of the final partons/particles and their fractional momenta xa,b is 

xa,b = —?= exp {±y). (22) 

At midrapidity at RHIC, only particle production with small p^ ^ 1 ^ 2 GeV will 
be sensitive to the shadowing region in parton densities. At the LHC, the region of 
transverse momenta will be much larger, see Fig. IHlfor the relevant kinematical regions 
which will be studied in proton-nucleus collisions at RHIC and the LHC for different 
observables. In ^2| extensive studies of the impact of parton densities on particle 
production with a large scale (e.g. transverse momentum or mass) at the LHC can be 
found. 

Most of the discussions at RHIC are done in terms of the so-called nuclear 
modification factors: 

D _ produced in A- B Na-b ,^^. 

expected in A- B N^oU X A^nucleon-nucleon ' 

Here Na-b is the number of produced particles of a given type in some kinematical 
region, in A-B collisions. Ncou is the number of binary incoherent nucleon-nucleon 
collisions with which particle production is expected to scale in collinear factorization, 
so Ra-b — * 1 for a large scale of the produced particle. 

In Fig. El the nuclear modification factor for neutral pion production in d-Au 
collisions at RHIC for pseudorapidity ?] ~ is shown and compared with results from 
collinear factorization |lU6j . Nuclear effects in deuterium, as stated previously, are 
usually neglected and, in any case, very small. From this Figure it becomes clear that at 
midrapidities only the region of small transverse momentum [px ^ 1 GeV) is sensitive 
to the amount of shadowing in parton densities. Above p^ ~ 3 GeV, one can see a 
(small) enhancement of particle production in nuclear targets - a phenomenon known 
as the Cronin effect for more than 20 years |135j . 

A large quark mass also provides a scale which gives some justification to the use 
of collinear factorization. Thus heavy flavour production can be employed to constrain 
parton densities inside nuclei |1361 1137] and their impact parameter dependence |138j . 
At RHIC, direct measurements of open heavy flavour production are still at the 
beginning jl39j . Most of the existing data are obtained from electrons from semi- 
leptonic decays |14()^ I141J which are weakly correlated with the parent heavy flavour. 
On the other hand, there exist studies of charmonium production. In Fig. El results 
for J/ip production in d-Au collisions at RHIC |142j for different pseudorapidities are 
shown and compared |143j to models for partons densities, using collinear factorization. 
Such models imply for forward rapidities, see (|^. an additional suppression (i.e. gluon 
shadowing as heavy flavour production is mostly sensitive to the gluon channel) on 
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Figure 16. Nuclear modification factor for neutral pions at 77 ~ in d-Au collisions 
at ^/s^ — 200 GeV. The curves are the results from |lU6j and the experimental data 
from [231 . [Figure taken from fTU^ .] 



top of the nuclear absorption of quarkonium in nuclei |144| I145J . From Fig. IT7I right, 
models which do not consider a very strong gluon shadowing jMl I103J are apparently 
favoured over those which consider a large shadowing |46j, see also the discussion at the 
end of Section IHl But the large error bars in the experimental data and the remaining 
uncertainty in the amount of nuclear absorption and its behaviour with rapidity |63| l64j. 
make it difficult to extract any definite conclusion. 

A striking finding in d-Au collisions at RHIC has been the change from 
an enhancement of the nuclear modification factor for light particles at central 
pseudorapidities (the Cronin effect, see Fig. ITI)jl to a suppression, at all measured 
transverse momenta, at forward pseudorapidities [TH IT^ IT^ I147J fFig. ITHjl ft. see the 
review |148j . In Glauber-like models as the MV model, the Cronin effect results 
naturally from the transverse momentum broadening due to multiple scattering. But 

ft The decrease of the nuclear modification factors for light and heavy hadrons with increasing rapidity 
has been a well known phenomenon for many years, but the transverse momentum structure of the 
suppression far from midrapidity has only recently been measured. 
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Figure 17. Plot on the left: rapidity distribution of the i/ip cross section. Plot on the 
right: nuclear modification factor for J/i/; |142| in d-Au collisions at ^s„„ = 200 GeV. 
The curves in this plot labeled Vogt |143| correspond to different nuclear absorption 
cross sections and different models for nuclear parton densities EKS O I103| and 
FGS 0ni) while that labeled Kopeliovich corresponds to [SJ- [Figure taken from |142| .] 



it has been taken as a great success of saturation physics that non-hnear small- a; 
evolution was able to predict the suppression at all transverse momentum from initial 
conditions which contained the Cronin effect [^ I149t \\bO\ 1151] . Studies within collinear 
factorization |152j indicate that shadowing has to be stronger than usually considered in 
DGLAP or Gribov inelastic shadowing |153j approaches, in order to justify these forward 
data, and that other effects are at work. In any case the clear conclusion can be drawn 
that such data require a large amount of shadowing, even more considering that the 
average values of x probed in the nucleus are not very small, (x) ~ 0.01 (assuming 2 — >■ 2 
processes) |152j j:. Other observables like Drell-Yan production pT j I154[ I155| I156| I157J . 
and midrapidity production of high-pj- hadrons in proton-nucleus collisions at the LHC, 
should further clarify this issue. 

d-Au collisions at RHIC have been crucial to establish the baseline on top of which 
final state effects due to a new state of deconfined matter, the Quark Gluon Plasma, 
can be searched for and eventually characterized |158j . Due to this possibility, the 
determination of nuclear shadowing in nucleus-nucleus collisions is far more complex 
than in proton-nucleus, as many other effects can be at work. In any case, the influence 



X This value of x lies at the upper border of applicability of saturation ideas, as commented in 
Subsection 12.31 Thus it is not clear that they can be employed to explain this phenomenon. 



Nuclear shadowing 



27 



oi, 



P 1 



11 = 



2 



^: 



^1 



• 0-20%/60-80% 
O 30-50%/60-80% 



11=1 



3>: 



^ 



2 



4- 



T| = 2.2 



O-O" 






_-c>o-0-(>- 



T1 = 3.2 



D-" 



.oo" 



.oo-O- 



12 3 4 5 12 3 4 

Pt [GeV/c] p^ [GeV/c] 



12 3 4 5 

Pt [GeV/c] 



12 3 4 5 

Pt [GeV/c] 



Figure 18. Ratio of nuclear modification factors for different centrality classes, for 
charged (two leftmost plots) and negative (two rightmost plots) particles at different 
pseudorapidities for central/peripheral (filled points) and semicentral/peripheral (open 
points) d-Au collisions at JT^ = 200 GeV |^. [Figure taken from TT.] 



of shadowing on bulk particle production will be large, resulting in a strong reduction of 
multiplicities in Au-Au collisions at RHIC and Pb-Pb collisions at the LHC, see e.g. the 
reviews |159tll6UJ . or |13mil3ltll61j . For example, in the approach |1H] based on Gribov 
inelastic shadowing and assuming some kind of factorization based on the AGK rules, 
reduction factors in head-on heavy ion collisions ~ 2 at RHIC and ~ 4 at the LHC can 
be expected for particle production at mid-rapidity, compared to p-p collisions scaled by 
the number of binary collisions. Other approaches [HIj based on saturation ideas which 
directly relate the shadowing measured in lepton-nucleus collisions to multiplicities in 
proton-nucleus and nucleus-nucleus collisions, also predict large suppressions and result 
in agreement with experimental data, see Fig. ^| Nevertheless it must be stressed that 
the uncertainties on the separation of effects in the nuclear wave functions from other 
effects due to the collision are large, as there is no well established factorization scheme 
presently available for nucleus-nucleus collisions. 



6. Conclusions 



The phenomenon of shadowing is of large importance from a theoretical point of 
view: the behaviour of the nuclear wave function at high energies provides useful 
information for our understanding of QCD in such regime. It has also very strong 
practical implications: parton densities in nuclei are required to predict and understand 
particle production in collisions involving nuclei. In this article the phenomenon of 
nuclear shadowing has been introduced. I have discussed multiple scattering as the 
underlying physical mechanism illustrated through a simple example of two scatterings, 
and presented several models which make use of it. Then I have analyzed the DGLAP 
approach which does not address the origin of shadowing but only the evolution 
of nuclear parton densities through the DGLAP equations. Next I have shown a 
comparison of the different models for F2 and xg. Finally, I have reviewed the results 
in nuclear collisions for which shadowing is expected to play a large role, with special 
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Figure 19. Energy and centrality dependence (in terms of the number iVpart of 
nucleons participating in the collision) of the multiplicity of charged particles at 
midrapidity divided by A'part, in Au-Au collisions at several RHIC energies and for LHC 



energies. The results ^ of (2/Np^,t)[dN^^^ /drf\{r^ - 0) = 0.47[^ 



■10.288 



iVO, 



.089 
part 



are 



compared to PHOBOS data |162[I163J . Also shown in the lower panel are the p-p data 
|164lll6i3] . as quoted in 11631 . [Figure taken from ^91^.] 



emphasis on proton(deuteron)-nucleus collisions at RHIC. 

The uncertainties on gluon shadowing at small x < 10~^ are as large as a factor 
3, see Fig. HSl Multiple scattering approaches suggest a large amount of shadowing 
and tend to indicate that shadowing for gluons is stronger than that for quarks. But 
DGLAP evolution disfavours such situation in view of the existing data on Q^-evolution 
of the nuclear modification factors, at least for x ~ 0.01 where experimental data lie. In 
this x-region the uncertainties in multiple scattering approaches are large due to finite 
coherence length effects. Besides, at the existing values of Q^ for small x the validity of 
pure DGLAP evolution or the need of power-suppressed corrections, re-summed in the 
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totally coherent limit in some approaches like high-density QCD, is not yet clear. So 
more data at smaller x, eventually coming from a large energy lepton-nucleus collider 
like the EIC |21!, or from UPC [11311111, are needed. 

d-Au collisions at RHIC provide very useful information if the validity of some form 
of factorization is assumed. J/^ production data in the forward (deuteron) direction 
constrain, within coUinear factorization, shadowing for gluons, apparently favouring 
models with not very large shadowing. On the other hand, nuclear modification factors 
for charged or negative particles in the forward direction show a suppression for all 
transverse momenta first predicted in the framework of high-density QCD. They seem 
to indicate that shadowing has to be larger than in most available models. Both sets of 
experimental data truly imply more physical mechanisms than nuclear shadowing, and 
some sizable uncertainty still comes from the use of any kind of factorization. 

Nucleus-nucleus collisions at RHIC and the LHC will be strongly affected by 
shadowing. Thus, the quantification of shadowing effects will be crucial for the eventual 
characterization of a dense medium where many additional physical processes may be 
at work. In this respect, in the near future the LHC will open new possibilities, both 
with the heavy ion program and, mainly, with proton-nucleus runs pTj . 
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